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Abstract. Based on distributions of local Green's functions we present a stochastic approach 
to disordered systems. Specifically we address Anderson localisation and cluster effects in binary 
alloys. Taking Anderson localisation of Holstein polarons as an example we discuss how this 
stochastic approach can be used for the investigation of interacting disordered systems. 



1. Introduction 

In many cases the influence of crystal imperfections present in any real material can be neglected 
in favour of the interactions dominating e.g. electron transport or optical properties. In some 
cases however, and especially at low temperatures, interesting physical effects arise from the 
specific physical processes induced by scattering on crystal defects. Examples range from the 
quantum Hall effect - an intricate problem far beyond the scope of our study - to electrical 
transport in polaronic systems. One prominent example is Anderson's prediction that in a 
disordered material itinerant (extended) states can be turned into localised states which do not 
carry any current pQ. This transition from extended to localised states is a result of quantum 
interference arising from elastic scattering on impurities in a crystal. It leads, roughly speaking, 
to spatial confinement of an electron: The electron wave-function decays exponentially with 
distance. 

This peculiar, essentially quantum mechanical behaviour has motivated lots of research (see 
e.g. [2j for a review). As yet, however, many important questions in the field of localisation 
physics have not been finally settled, and localisation in interacting systems is still far from being 
understood. Progress on this topic requires new techniques which allow for a comprehensive - 
necessarily approximate, but reliable - description of disorder and interaction on a microscopic 
scale. One possible approach shall be discussed here. 

The focus of our study is on substitutionally disordered three-dimensional materials like doped 
semiconductors or alloys. We will not be concerned with amorphous materials like glasses which 
do not possess a crystal lattice (but see e.g. (HI El)- Then disorder primarily manifests through 
site-dependent local potentials ej. The motion of an electron in such a disordered crystal is 
described by the tight-binding Hamiltonian 

H = j2^4 c i- t J2 c h- w 

* (iJ) 



In a perfect crystal, i.e. ej = 0, the electron hopping t between nearest-neighbour sites 
gives rise to a band of width Wo, e.g. on a cubic lattice Wq = 12t, or Wq = AtVK on a Bethe 



lattice with connectivity K. In what follows, we fix the bandwidth Wq = 1, and measure energies 
in units of Wo. 

To model the stochastic character of disorder the ej are considered as random variables 
with a given distribution. We demand that they are independently identically distributed with 
common distribution P(ej). For such a stochastic Hamiltonian quantities like the (retarded) 
Green's function Gy(cj) = lim r? _ >0 + (0| c « [ u + — H]~ 1 c j|0) are random variables in their own 
right. It is thus reasonable to ask for their distributions. This is especially true for the local 
density of states (LDOS) 

Pi(u) = -Im Gu(uj)/ir . (2) 

Its distribution P[pi(u)] gives the statistics of the LDOS in a disordered system, where 
translational symmetry is broken and Pi(to) varies with i. Note that P[pi{u))] does not depend on 
the site index i: Translational symmetry is restored on the level of distributions. It is important 
to realise that the site-dependence of the LDOS constitutes an eminent aspect of a disordered 
system. For an extended state, for example, Pi(u>) has a finite value on most lattice sites. For a 
localised state, in contrast, Pi(to) has exponentially small values on lattice sites outside a finite 
region, corresponding to the decay of the wave- function. P[pt(u)] will thus be different for 
extended and localised states. Obviously the disorder averaged density of states (DOS) 

oo 

PaveM = (pi(oj)) = J Pi P[pi((jj)] dpi , (3) 


merely counting the number of states at a given energy u, cannot account for this signature. 
Only for the ordered case P{ei) = <5(ej), when translational symmetry implies that pi does not 
depend on i, the distribution is entirely determined by the DOS, P[pi(u)] = 5[pi — p ave (u>)]. As 
we will later see a similar characterisation does not even hold approximatively in a disordered 
system. The full distribution can be extremely broad for states in which the electron strongly 
scatters at impurities, catching the resulting deviations in pi. 

2. Local distribution approach 

Our goal is to obtain a calculational scheme for the distribution of the LDOS for the models given 
by Eq. This scheme must account for the correlations that make up Anderson localisation. 
At the same time it should be extendable to incorporate interactions. Naturally this requires 
to employ approximations. It is important to guarantee that within these approximations the 
mutual influence of interaction and disorder at the microscopic scale is still represented (which 
rules out e.g. the coherent potential approximation (CPA), see below). A definite candidate to 
meet these demands is what we call the local distribution (LD) approach. It provides a self- 
consistent scheme for the distribution P[Ga(u>)] of the local Green's function Gu(uj). We first 
describe this approach and its application to disordered systems, and later will combine it with 
a treatment of interaction. 

2.1. LD approach on a Bethe lattice 

It is straightforward to construct the LD approach following the work of Abou-Chacra, Anderson 
and Thouless 0. We briefly repeat their derivation which is carried out on a Bethe lattice (see 
Fig. The basic observation is the decomposition of the local Green's function Gu(lv), 
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Figure 1. Sketch of the (half- infinite) Bethe 
lattice with coordination number K = 2. When 
site i is removed from the lattice the sites j are 
in the same geometrical situation as site i before, 
and no path connects them. 



The sum runs over the K neighbours j of i, and the Green's functions Gjj (co) have to be 

calculated for the lattice after site i is removed. For the ordered case we have G§ = Gu{u) (cf. 
Fig^), and Eq. Q) is a quadratic equation for the Green's function. Especially we can read off 
the semi-circular DOS p(u) = 4/(ttW 2 ) - 4w 2 of the Bethe lattice 

Now the stochastic viewpoint comes into play. For a disordered system, like all Green's 
functions are stochastic quantities, i.e. random variables, and Eq. (JIJ determines the random 

variable Gh(uj) in terms of K + 1 others, namely and G^{uj). On the Bethe lattice the 

Green's functions G^(uj) correspond to the same geometrical situation as Gu(lu), hence they 
are identically distributed. Furthermore no path connects the sites j = 1, . . . , K if i is removed. 
This implies that the Gjj (to) are independently distributed. Therefore Eq. Q can be regarded 
as a self-consistency equation which determines the distribution of one random variable, the local 
Green's function Gu(ll>), or equivalently G^(uj). The distribution P(e«) enters this equation as 
a parameter. 

Such a stochastic self-consistency equation can be numerically solved by a Monte-Carlo 
procedure (Gibbs-sampling). To this end we represent the random variable Gh{oj) through a 
sample of N entries (typically N = 10 4 . . . 10 7 ). Each entry of the sample is repeatedly updated 

by a new value from Eq. ([4"]). with K values for G^J(lo) randomly drawn from the sample and 
a randomly chosen e,. The distribution of Gu(uj) is finally constructed like a histogram by 
counting the number of entries in the sample with a specific value. 

2.2. LD approach on arbitrary lattices 

To put the LD approach in a more general context we will show how to construct it as 
an approximative scheme on others than the Bethe lattice, thereby making contact to the 
CPA (concerning CPA see e.g. 0). On an arbitrary lattice the Green's function can still be 
decomposed as 



In contrast to the Bethe lattice the Gf- (u) which appear in this decomposition do not correspond 
to the same geometrical situation as Gu(uS). If we iterated Eq. Q we would thus end up 
with an infinite hierarchy of equations containing many different Green's functions. As an 
immediate consequence the LD approach, relying on local Green's functions, will provide only 
an approximative scheme on general lattices. 

A simple approximation is suggested by introducing 
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which is the part of the K hopping contributions in Eq. Q where the electron leaves site i 
through j. Specifically we assume that, for each j, H^\uj) is determined by the local Green's 
function G 



alone. Moreover, the K Green's functions G 



are taken to be independently 
distributed. 

To proceed further we write the local Green's function with a self-energy like in the CPA 

Gu(u) = G°(lo- EiH) , (7) 

where G°(uj) is the 'bare' propagator for the ordered system. In contrast to the CPA the 
self-energy is now site-dependent. Thinking in terms of an effective medium in which 

site i is embedded this medium is characterised by the distribution of self-energies 
The electron at site i 'sees' the effective medium through the K Green's functions Gjj(co). 

Being part of the effective medium the hybridisation H (co) associated with Gjj(uj) must fulfil 
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Inverting this equation yields 



Hf{u) = ^(u-^-[G 3 Mr 1 



(8) 



The Green's function Gu(lo) is then given by K other Green's function Gjj(ui) of the same 
type: 
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(9) 



Clearly, if all = 0, this equation reduces to £j = 0, that is Gu(lo) = G% (u). 

At last, we have a complete set of equations, which with the same Monte-Carlo-scheme as 
before can be solved to determine the distribution of Gu(lu). Note that the approximations we 
made are guaranteed to respect causality, in the sense that lmGu(u>), ImEj(w) < 0. Eqs. |Zj) 
and © are immediately seen to preserve this important property. For Eq. Q, where this 
property is not obvious, we can rely on a result from the analyticity proof of the CPA 8 , which 

states that Im |s + [G°(u - £)] _1 } > for ImE < 0. 

We like to point out that the construction done here is slightly ambiguous. A different 
approximation in Eq. (jHJ would still give a self-consistent scheme for distributions of local Green's 
functions Gh(uj). Anyhow, our scheme incorporates the CPA and reproduces the exact scheme 
for the Bethe lattice in a natural way. The CPA results from Eq. @ for K = oo by the central 
limit theorem. Then, the sums over Green's functions and self-energies are replaced by averages, 
yielding the CPA equation for the averaged local Green's function 
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(10) 



In principle, K is in our scheme no free parameter but the number of neighbours to a lattice site. 
If we increase K the bare propagator G°(uj) changes, and therefore we shall scale t as t oc l/y/~K 
like in the limit of high dimension Then, for K — > oo, our scheme still reduces to the CPA, 
which is known to become exact for lattices with infinite connectivity JHJ- This implies that 
the approximations we made are good for high-dimensional lattices. Actually they turn out 
to be good for d > 3. It is commonly known that CPA is the best single-site approximation 
for disordered systems. The local Green's function is obtained on average from an 'averaged' 
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Figure 2. Change of the 
probability distribution P[pi(u>)] 
of the LDOS Pi(u>) for the 
Anderson model, in the band 
centre u = 0, with increasing 
7. Note that the inset has a 
logarithmic /Oj-axes. 



effective medium. The LD approach solves the local problem exactly, while the hybridisation 
is calculated from a 'fluctuating' effective medium. So to say, the LD approach is the best 
single-site approximation if fluctuations are taken into account. 

On the Bethe lattice the Green's function fulfils G^(u) = [u — ^KG^lu)] (this is just 
Eq. (@J) for t{ = 0), whereby Eq. (jJJ) reads 

£4(0/) - [Guiu)}- 1 — uj fKGniuj) (11) 

Then Eq. © reduces to Eq. Q for the Bethe lattice. The key point is of course that, owing to 

(i) (i) 

the specific Bethe lattice geometry, the non-diagonal contributions G^ are zero, and the Hj 
are indeed independently distributed. Then our approximation scheme becomes exact: On the 
Bethe lattice the best single-site approximation, with fluctuations, is an exact theory. Note that 
for the Bethe lattice K is a free parameter: The 'free' DOS in absence of disorder is, for fixed 
Wo, independent of K, but choosing K fixes the Bethe lattice used. 

As the transfer of the LD approach to other lattices shows, going to the Bethe lattice 
basically changes the 'free' DOS and introduces approximations to cubic lattices by neglecting 
correlations. The essential properties of disordered systems, e.g. localisation, will however be 
well described. For the time being we work on a K = 2 Bethe lattice. 

3. Disordered crystals 

Below we will discuss two models, the Anderson model with a 'continuous' probability 
distribution, and the binary alloy model with a 'discrete' one. While the former is dominated by 
the Anderson transition from extended to localised states, the latter is dominated by multiple 
electron scattering on clusters of atoms, visible e.g. in the fragmentation of the DOS (see below). 

3.1. Anderson model 

In the Anderson model, with disorder strength 7 > 0, the e« have a box distribution 

P(ei) = (l/ 7 ) 6(7/2 -N). (12) 

Fig. El shows the change of the distribution P[pi(uj)] of the LDOS with increasing disorder, 
as obtained from our Monte-Carlo scheme. We find the behaviour indicated in the introduction. 
For small disorder the distribution is rather Gaussian, but becomes strongly asymmetric 
and extremely broad as 7 increases. The (averaged) DOS p &ve {u) is not indicative of these 
distributions. If 7 exceeds a critical value 7 c (w) [7 c (w = 0) ~ 3.0] the distribution becomes 
singular while p ave (uj) is still finite. At this point the transition from extended to localised 
states takes place. To determine 7 c (w) numerically we have to study the scaling of appropriate 




Figure 3. The figure displays, 
for localised states (7 = 1.5, 
uj = 0.9), the change of P{pi\ 
with decreasing 77. The inset 
shows that p ave is constant for 
77 — >• 0, but a typical moment, 
the geometrically averaged DOS 
p geo = exp(lnpj), goes to 
zero. For extended states we 
get distributions like in Fig. [21 
independent of (a small) 77. 



quantities, like for any phase transition. Here we employ the different behaviour of P[pi(uj)\ 
when the imaginary part 77 in the energy argument of Ga{u + irj) goes to zero: For extended 
(localised) states the distribution is stable (instable) for 77 — > (see Fig- EI) - This behaviour is 
again not caught by p ave (w), which is finite for 77 — ► for both extended and localised states. But 
suitable moments of the distribution, like the geometric average, drop to zero just for localised 
states (see inset in Fig. |3J). With this criterion which numerically tests whether the distribution 
is singular or not we can very precisely distinguish extended and localised states Fig. 0] 

shows the phase diagram of the Anderson model on the Bethe lattice obtained by means of 
the limit 77 — ^ 0. It displays two important features of localisation. First, the existence of 
a critical disorder j c (lo = 0) above which all states are localised due to the strong impurity 
scattering. Second, before complete localisation occurs, states towards the band centre are 
extended while states towards the band edges are localised, these states being separated by 
the so-called mobility edges. Furthermore the mobility edge trajectory reveals that for small 
disorder the electron tunnels between impurities, giving rise to extended states outside the band 
[—Wo/2, Wo/ 2] of the ordered system. 

Let us finally contrast the LD result for /9 ave (u;) (Fig. EJ) with the CPA. The CPA reproduces 
the DOS very good inside the band - only here P\pi(uj)\ is a Gaussian - but fails closer to the 
band edges where it misses the (Lifshitz) tails in the DOS. These tails result from states at the 
band edges which come along with repeated scattering of the electron on clusters of impurities 
with a large (or small) potential e». Since multi-scattering is not accounted for in the CPA these 
states cannot be described. 

3.2. Binary alloy 

Multi-scattering becomes rather important for 'discrete' distributions like 

P(ei) = c A 5(€i + A/2) + (1 - c A )5(€i - A/2) , (13) 

describing a two-component binary alloy made of 'A- atoms' with potential —A/2 and 'B-atoms' 
with potential +A/2. Thus we have two parameters, the concentration ca of the A-species and 
the energy separation A. 

Exemplarily we consider the case ca = 0.1 when the A-atoms form the minority species in a 
bulk of B-atoms, and A = 2.0 when the A-atoms are energetically well-separated from the band 
of the B-bulk. Fig.EJshows the corresponding DOS, contrasting the LD and CPA results. The B- 
band is rather smooth and almost reproduced by the CPA, which nevertheless misses important 
features. In the energy range of the A-atoms we do not find a band but a strongly fragmented 
set of peaks. Each of these peaks can be attributed to a specific cluster of A-atoms, as indicated 
in the figure. Due to the low concentration ca of A-atoms, and the large energy separation 




Figure 4. Left panel: Mobility edge trajectory, i.e. the phase diagram, of the Anderson 
model, obtained in the limit rj — ► 0. The dashed lines indicate the band edges ±(W/2 + 7/2). 
Right panel: DOS for the Anderson model at 7 = 1.5 for oj < 0, calculated with CPA and LD 
approach. Both the DOS and the phase diagram are symmetric under to — > —oj. 

A, states at these clusters do not hybridise and form a continuous band, but are essentially 
confined to the clusters and contribute discrete peaks to the DOS. The CPA is by construction 
unable to identify those cluster states and therefore entirely misses band fragmentation. With 
increasing A these signatures become more prominent, while CPA suggests that the bands, being 
well separated, less influence each other (see right panel in Fig. [5J). The CPA is of course able 
to depict gross features as band splitting, but replaces the complex spectrum by semi-circular 
bands with weight ca and 1 — ca- Then, e.g., the minority A-band has spectral weight outside 
the CPA band: Replacing this band by a semi-circular one while conserving the overall spectral 
weight ca results in a clearly reduced CPA bandwidth. 

We think that the binary alloy is a particular instructive example how, using the distribution 
of Gh{uj), one manages to account for multi-scattering events which are absent in a CPA 
description relying on averaged values. With a proper treatment of multi-scattering like in 
the LD approach both quantum interference leading to Anderson localisation and formation of 
cluster states is correctly described. 

4. Electron-Phonon Coupling 

Anderson localisation arises from elastic electron-impurity scattering. Coherence which is 
maintained during elastic scattering is the important precondition for localisation. To 
understand which role localisation plays in a real solid we must therefore understand how it 
is affected by inelastic scattering, which naturally arises from the interaction of the electron 
with lattice vibrations, i.e. phonons. A model to study the competition between localisation - 
related to coherence- and electron-phonon (EP) interaction -related to incoherence- is provided 
by the Anderson-Holstein Hamiltonian 

H = e 4 c * - t Yl c h - v^J2( b t + b ^ c l ci + ^Y1 b l b i > ( 14 ) 

which adds a Holstein-type of interaction to the Anderson model, the being distributed 
according to Eq. (|T2|) . Here denote bosonic operators describing dispersionless optical 
phonons with frequency u>q which are locally coupled to the electron density cfcj. We introduce 
the dimensionless phonon frequency Cjq = ujq/Wq, and a coupling constant A = 2e p /Wq. 
Moreover we consider one electron at T = 0. 



Q. 




Figure 5. Left panel: DOS for the binary alloy model at A = 2.0, ca = 0.1, contrasting the 
CPA and LD results. Arrows mark different peaks in the DOS which arise from the clusters 
indicated (e.g. A-BB indicates the contribution from an A-atom which is surrounded by two 
B-atoms.) We used in the figure an arbitrary r\ = 10 -3 to broaden the peaks. Right panel: Kind 
of a phase diagram for the binary alloy model at ca = 0.1, showing the DOS for various A. The 
dashed curves show the CPA band edges, and the dotted lines mark u> = ±A/2 ± Wo/2 (figures 
taken from [T2] ) 



Besides introducing incoherent motion, EP coupling in this model can lead to the formation of 
a new quasiparticle. Sufficiently strong EP coupling binds the electron to the lattice deformation 
at the electron's site, forming a new compound entity, the polaron (see e.g. |13j). A polaron 
is strongly mass-enhanced and therefore very much affected by disorder. However, a polaron 
is not just a heavy electron. Due to inelastic scattering and retardation of the EP interaction 
in most cases the internal structure of the polaron will play a crucial role, leading to different 
localisation properties. Moreover if inelastic scattering dominates, leading to incoherent motion 
of the polaron, localisation will be suppressed. Then a polaron can be less easy to localise than 
the free electron, although its mass is still increased. 



4-1- Extending the LD approach to interacting systems 

Like in the standard Green's functions formalism interaction is incorporated in the LD approach 
by an interaction self-energy H^(u) The 'disorder' self-energy Sj, as introduced in Eq. (JJJ), 
is local. The interaction self-energy will be approximated to be local as well. The Green's 
function Ga{u) is then obtained as 
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(15) 

replacing Eq. ©. Like for the disorder part one should demand that the interaction self- 
energy is obtained as the 'best' local self-energy which amounts to use the dynamical mean field 
theory (DMFT) (for a review see |16j). Within DMFT the interaction self-energy is a functional 
S[ (uj) = S[ [-Fj(w)] of the local propagator 
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which does not contain £^(u>). With the DMFT self-energy inserted into the equations still the 
Monte-Carlo scheme applies. Since the interaction couples different energies each entry of the 
sample now represents a local Green's function Gu(uj) on a set of u (e.g. Matsubara frequencies). 
Note that during each update of an entry one has to evaluate £f[-Fj(u>)]. Without disorder, for 
€i = 0, this scheme reduces to standard DMFT. 

The functional dependence £^[^(0;)] is not explicitly known but requires, despite the 
approximations made, the solution of a quantum-mechanical many-particle problem. Solving 
this model constitutes the main part of a DMFT implementation. Concerning disordered 
interacting system the evolved calculation of X^fi^u;)] poses a severe problem. So far system 
with a finite charge carrier density like a disordered Hubbard model could only be addressed 
in limiting cases |15j . For a single electron in the Anderson-Holstein model we are fortunate to 
explicitly know the functional [Fi(co)], given as an infinite continued fraction for T = |17U18j : 

= 1£PUJ ° or- • (17) 
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The p-th level of this continued fraction corresponds to the emission of p virtual phonons, shifting 
the energy argument of Fi by pco®. Originally this expression has been obtained as an extension of 
the CPA to describe dynamical interaction effects (concerning dynamical CPA, see JZ])- I* 1 the 
spirit of the CPA the interaction can be mapped onto an effective Hamiltonian which is related 
to the original problem via a consistency condition on the Green's function |19j . This mapping 
can be combined with the LD construction on the Bethe lattice whereby the consequences of 
the approximations concerning interaction and disorder become very explicit It turns out 
that the interaction induced correlations between local Green's functions are almost treated on 
the same level of correctness as the disorder induced correlations. In particular one finds that 
the extended LD approach comprises cooperative effects which arise from the mutual interplay 
between interaction and disorder. This is why the LD approach excels approaches which try 
to replace the interacting disordered system by an effective non-interacting disordered system. 
The interaction is there mimicked by effective transfer integrals or local potentials, but the 
feedback of disorder on interaction is not accounted for. The same objection applies if the 
interacting disordered system is replaced by an effective interacting system, mimicking disorder 
by effective coupling constants. The price to be paid on using the LD approach is the rather 
involved numerical implementation combining a repeated solution of the DMFT problem with 
a Monte-Carlo scheme. 



4-2. Polaron localisation 

Within the extended scheme just described we can study how the electron is affected by disorder 
and EP interaction. We will only touch on this question and focus on two significant cases where 
we compare parts of a phase diagram for polaron localisation to that of the Anderson model. 
The reader should be aware that we encounter a very complicated physical situation. Already 
without disorder the Holstein EP interaction displays rich, and very distinct, physics in different 
parameter regimes. With disorder, cooperative effects like the formation of polaron defect states, 
come into play. These effects do hardly fit into a 'universal' phase diagram of polaron localisation 
but demand a more concrete description, depending on the specific case studied. A more detailed 
discussion of this and related issues can be found in |2Uj . 

In Fig. the DOS for two choices of polaron parameters as obtained by DMFT is shown. 
We used the 77 — * — criterion to obtain the mobility edges for the lowest polaron sub- 
band in each of the two cases. The comparison to the mobility edges of the Anderson model 




CO CO 

Figure 6. DOS p{oS) without disorder (7 = 0), obtained within DMFT. ImS(w) is displayed 
downwards. Left panel: The strong-coupling, large phonon-frequency case A = 9.0, u)n = 0.5625 
can be understood from the atomic limit (independent boson model), treating t as a perturbation. 
The DOS shows narrow bands which are separated by ujq. The bandwidth of the lowest band 
(magnified in the bottom row) is W/Wq = 3.45 x 10 -4 , i.e. the polaron being extremely heavy. 
Right panel: For weak coupling (u)n = 0.05, A = 0.25, upper row) EP interaction essentially 
leads to inelastic motion of the electron (ImS(w) < 0). For stronger coupling (A = 1.0, same 
u)n = 0.05) a polaron has formed. The lowest polaron band {W/Wq = 8.123x 10~ 3 ) is asymmetric 
revealing the different character of the polaron to the lower and upper band edge. 
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Figure 7. In addition to the 
phase diagram of the Anderson 
model the mobility edges for the 
two polaron bands from Fig. 
are shown. The mobility edges 
are drawn relative to the sub- 
band position, and disorder is 
measured in units of the sub- 
band width W |2Q|. 



(Fig. EJ) demonstrates that only in special cases the localisation properties of the polaron can be 
understood as the mere result of the mass renormalisation of the quasiparticle. 

For strong-coupling and large phonon-frequency A = 9.0, cDn = 0.5625 (circles in Fig. [7J), 
the polaron in the lowest band is basically a heavy particle which moves fully coherent but 
within an extremely narrow band. It localises like a bare electron with a rescaled bandwidth, 
and the internal structure of the polaron does not play a role. Due to the strong coupling the 
relevant energy scale for localisation has however changed by four orders of magnitude, making 
the polaron extremely susceptible to disorder. 

For intermediate coupling and small phonon-frequency A = 1.0, ojq = 0.05 (squares in Fig.[7J), 
the polaron motion is, already without disorder, different at the lower and upper band edge. 
At the lower band edge the polaron is rather mobile, while at the upper band edge it tends 
to be immobile. Concerning the localisation properties the polaron at the lower band edge is 
thus almost unaffected by disorder and does not readily localise. At the upper band edge, in 
contrast, it is easily localised. Note that the lowest polaron band is again coherent: localisation 
is not just affected by incoherent motion but intricately depends on the internal structure of 
the polaron. Be also aware that disorder is measured in units of the renormalised band width: 



the polaron at the lower band edge is more difficult to localise than a bare electron only with 
respect to the relevant energy scale W, i.e. the polaronic bandwidth. 

In Fig. 0we have only shown data for moderate disorder, i.e. below j c (u = 0) when, for 
the bare electron, all states become localised. In the first (strong-coupling) case we can indeed 
localise all states in the lowest polaron band and obtain the respective 7 c (w = 0) which has the 
value of the bare electron renormalised by W/Wq. Here disorder induced mixing with excited 
polaron states is prevented by the large gap to the next band (W <C uo). In the second (adiabatic 
intermediate-coupling) case, however, with increasing disorder the lowest band merges with the 
next polaron band before it is completely localised. Then the relevant energy scale for 7 changes 
by one order of magnitude to the joint bandwidths of the two lowest bands. Hence an equivalent 
to 7c(w = 0) does not exist. Apparently we can only draw parts of a phase diagram for polaron 
localisation. Again, EP interaction changes the localisation properties in a more complicated 
way than thinking only in terms of mass renormalisation would suggest. 

5. Summary 

In disordered systems Green's functions emerge as stochastic, random quantities. It is reasonable 
to take this stochastic character serious and to incorporate distributions of those quantities in 
a description of disorder. Following |Sj the proposed LD approach is based on distributions of 
the local Green's function. Being exact on a Bethe lattice, approximate on general lattices, it 
accounts for quantum interference leading to Anderson localisation, and multiple scattering on 
clusters like in a binary alloy. 

The LD approach is an extension of the CPA including fluctuations. Conversely, it contains 
the CPA in the limiting case when fluctuations are replaced by averages. We can then understand 
that the CPA results are good if - and only if - the distribution of the LDOS is characterised 
by its average. This condition does not mean that the disorder is weak: Even for small 7 in the 
Anderson model states at the band edges are localised, and beyond the scope of the CPA. A 
similar restriction exists for the minority band in the binary alloy model. 

The great challenge is however not disordered, but interacting disordered systems. Here 
the LD approach can provide a description of the microscopic interplay between interaction 
and impurity scattering. Locally, both interaction and impurity scattering are treated in the 
best approximation possible: DMFT is the best single-site approximation to interaction; and 
impurity scattering is exactly treated (on general lattices at least locally exact). 

The LD approach depends on the accuracy of the DMFT self-energy. A deep understanding 
of this part of the problem is important. How the involved numerical methods (quantum 
Monte-Carlo, exact diagonalisation, density matrix renormalisation group, to mention a few) 
or approximative schemes (e.g. iterative perturbation theory (IPT)) used for the calculation of 
the DMFT self-energy ^Sj cope with the considerable fluctuations in the local propagator Fi(tu) 
induced by disorder is by no means clear (e.g. the perturbative expansion in the IPT is best for 
weakly varying Fi(oj)). Since we largely understand the LD approach's application to disorder, 
how to guarantee for the quality of these methods is still the biggest problem - the computational 
costs for their implementation in the Monte-Carlo scheme notwithstanding. Our results for the 
localisation of a single Holstein polaron, where these objections are absent, demonstrate that the 
LD approach by itself allows for a comprehensive description and detailed study of interacting 
disordered systems. Physically we learned how intricate and rich the interplay of interaction 
and localisation can be. Future work will certainly support and extend this observation. 
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